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Abstract. In this article, we will briefly review the recent progress on hydrodynamic modeling and the 
extraction of the quark-gluon plasma (QGP) specific shear viscosity with an emphasis on results obtained 
from the hybrid model VISHNU that couples viscous hydrodynamics for the macroscopic expansion of the 
QGP to the hadron cascade model for the microscopic evolution of the late hadronic stage. 
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PACS. 25.75.-q Relativistic heavy-ion collisions 
flow - 24.10. Nz Hydrodynamic models 



12.38.Mh Quark-gluon plasma - 25.75.Ld Collective 



1 Introduction 

Heavy ion experiments at the Relativistic Heavy- Ion Col- 
lider (RHIC) at Brookhaven National Laboratory and the 
Large Hadron Collider (LHC) at CERN, have accumu- 
lated strong evidences for the creation of the quark gluon 
plasma (QGP) [l][2l[3]. The observation of strong collective 
flow and the successful descriptions from hydrodynamics 
demonstrate that the QGP is strongly coupled and be- 
haves like an almost perfect liquid with a very small spe- 
cific shear viscosity [HE]. Using gauge/gravity (AdS/CFT) 
correspondence, Kovtun, Son and Starinets showed that 
there exists a lowest limit for the shear viscosity to en- 
tropy density ratio rj/s = 1/Att (called as KSS bound) for 
a large class of strong-coupled quantum field systems (not 
including QCD) [5]. This raises the question how close to 
this limit is the specific shear viscosity of the QGP created 
at RHIC and the LHC. 

Using weakly coupled QCD, one can calculate the QGP 
shear viscosity at very high temperatures from kinetic the- 
ory or from Kubo formula [7]. However, it is difficult to do 
a first principle calculation for the shear viscosity of the 
strongly coupled QGP created at RHIC and the LHC. It 
thus desirable to extract it from experimental data. Pre- 
vious studies revealed that the anisotropic flow generated 
in relativistic heavy ion collisions is highly sensitive to the 
shear viscosity due to the rapid expansion of the QGP 
fireball, which leads to large viscous corrections from the 
shear velocity tensor 0. Explicit viscous hydrodynamic 
simulations revealed that even the small specific shear vis- 
cosities at the KSS bound leads to significant suppression 
of elliptic and triangular flows [ilinEIllIlliaillllSEn] . 
In principle, this allows for an extraction of the specific 
QGP viscosity rj/s by tuning rj/s in viscous hydrodynamic 
calculations and fitting the results to the sensitive exper- 
imental observables. In practice, this procedure requires 
sophisticated theoretical modeling of the heavy-ion colli- 



sions, such as the initial conditions that include the fluc- 
tuation effects, an equation of state (EOS) that properly 
describes the speed of sound and the QCD phase transi- 
tion, non-equilibrium kinetics and chemical composition 
in the late hadronic stage, etc. [T7] . 

For a realistic description of the late hadronic stage, we 
developed the VISHNU hybrid model that couples viscous 
hydrodynamics with a hadron cascade model which mi- 
croscopically simulates the hadronic rescattering and the 
chemical and thermal freeze-out of varies hadron species 
through solving the Boltzmann equation with flavor de- 
pendent hadronic cross sections |18j . In this article, we 
briefly review the recent progress in hydrodynamic mod- 
eling and the extraction of the QGP viscosity from elliptic 
flow data with a special emphasis on the results obtained 
from the VISHNU hybrid model [1^112011211122] . 



2 Viscous hydrodynamics and the hybrid 
approach 

2.1 Viscous hydrodynamics 

Viscous hydrodynamics is a macroscopic tool to describe 
the expansion of the QGP and subsequent hadronic mat- 
ter. In this section, we briefly review the Israel-Stewart (I- 
S) viscous hydrodynamics in 2+1 dimension with the as- 
sumption of longitudinal boost invariance, which is solved 
by the VISH2+1 code developed at the Ohio State Univer- 
sity around 2007 [TU]. For the work related to Ottinger 
and Grmela viscous hydrodynamics and recent progresses 
on 3-|-l-d viscous hydrodynamics (based on I-S formal- 
ism) without longitudinal boost invariance, please refer to 
Ref m] and Ref. [16l[23l[24| respectively. 

VISH2+1 solves the equations for energy-momentum 
conservation and the 2nd order Israel-Stewart viscous equa- 
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tions [251126] (For simplicity, the net baryon number and 
heat conductivity are assumed to be zero). With Bjorken 
approximation [27j . these equations can be conveniently 
written in the curvilear coordinates x"^ = {T,x,y^ris), 
where r = and 77,, = i In(f^) [IHlllH]: 

d„T™" = 0, T'"" ^ eu"'u" - pA""' + tt"", (1) 

-2" — ^^H^^j' 

D7T=- — (7T + Ce)-i7T^4fe'V (3) 
Tn 2 Tn \Qi J 

Here e is the local energy density, p is the local pres- 
sure, and m"* is the flow 4- velocity. U is bulk pressure and 
^mn jg i-jjg shear stress tensor. D = u^dm and V™ = A"^'di 
^^mn _ gmn_^m^n-^ ^-^^ -j-j^g ^^jjj^g spacial derivative 
in the local comoving frame, ct™" = i(V'"u"+V"u™) — 
i^mn^^^yfc-j ^j^g velocity shear tensor. 

The shear viscosity 77, bulk viscosity Q and the corre- 
sponding relaxation times and tjj are free parameters in 
viscous hydrodynamic calculations. The default settings in 
VISH2+1 are: 77/s = const., = iil/{sT) [10| and C/s = 
(We found that the bulk viscous effects are much smaller 
than the shear viscous effects due to the critical slowing 
down near phase transition [29]). For hydrodynamic cal- 
culations with temperature-dependent rj/ s{T), please refer 

to puin] 

The equation of state (EOS) is an additional input for 
hydrodynamic simulations. The default EOS s95p-PCE in 
VISH2+1 implements recent lattice results for the QGP 
phase and emphasizes the partially chemical equilibrium 
in the hadronic phase for temperatures below Tchom = 
165 MeV [3l[33]. 

The initial entropy density profile for VISH2+1 are pro- 
vided by two popular initial geometric models: Monte Carlo 
Glauber Model (MC-Glauber) and Monte Carlo KLN Model 
(MC-KLN) [33]. Due to the finite number of coUiding 
nucleons, the initial eccentricities (which are the driving 
forces for the elliptic flow) fluctuate from event to event 
for a specific central! ty bin. To account such fluctuating 
effects on average, we generate a large number of initial 
entropy density proflles from MC-Glauber or MC-KLN 
model, rotate each distribution either by aligning the par- 
ticipant plane or the reaction plane, and then average such 
rotated profiles to obtain one smoothed initial entropy 
density profile with participant plan eccentricity Spart or 
reaction plan eccentricity e^ec- For calculation efficiency, 
early VISH2+1 and current VISHNU perform calculations 
with such event-averaged initial conditions, which is called 
as single-shot simulations (With Spart that is significantly 
larger than e^ec in most central and most peripheral col- 
lisions, single-shot simulations with smoothed initial pro- 
files averaged in the participant plane partially accounts 
the fluctuation effects imprinted in the flow data measured 
in the participant plane, such as W2{2}). The fluctuating 
proflles can also be directly put into VISH2+1 and VISHNU 
in the event-by- event simulations, resulting in fluctuating 



hadron spectra and flow that varies from event to event, 
which are then averaged to compare with the experimental 
data. Currently, event-by-event simulations have not been 
implemented in VISHNU. For recent progress from event- 
by-event VISH2+1 simulations, one can refer to Ref. [151 

[sg. 

The decoupling temperature that defines the hydrody- 
namic freeze-out surface is generally set to T^cc = 100 — 
120MeV to allow for sufficient evolution time to build up 
the phcnomcnologically required radial flow which con- 
trols the slopes of the hadron spectra and their depen- 
dence on hadron masses |33j . 

2.2 VISHNU hybrid model 

Although the implementation of the EOS s95p-PCE prop- 
erly accounts for the chemical freeze out for various hadron 
species, the pure hydrodynamic approach ultimately fails 
in the late hadronic stage due to the dramatic increases 
of viscous corrections which invalidate the fluid dynami- 
cal approach that requires near equilibrium. For a more 
realistic description of the evolution and decoupling of 
the late hadronic stage, we developed the hybrid model 
VISHNU [18] by combining viscous hydrodynamics for the 
QGP fluid expansion with the hadron cascade model for 
the kinetic evolution of the hadronic resonance gas at a 
switching temperature Tsw 

The viscous hydrodynamics implemented in VISHNU 
is VISH2+1, which has been briefly described in Sec. 2.1. 
The hadronic cascade model used there is UrQMD (Ultra- 
rclativistic Quantum Molecular Dynamics model) |36j , which 
microscopically simulates the evolution of the hadron reso- 
nance gas through the coupled Boltzmann equations with 
flavor-dependent cross-sections. The connection between 
VISHNU and UrQMD is realized through a Monte-Carlo event 
generator called H20 which converts hydrodynamic out- 
put into particles proflles for further UrQMD propagation 
by sampling the Cooper- Frye phase-space distribution (in- 
cluding viscous corrections) on the decoupling surface |18) . 
The partially hadronic chemical equilibrium is naturally 
imprinted in UrQMD by simulating the dynamics of the 
hadronic gas with elastic, semi-elastic and inelastic colli- 
sions. By describing the hadronic rescattering and freeze- 
out procedure microscopically, VISHNU improves purely 
hydrodynamic models and eliminates the additional ad- 
justable parameters required for the transport and freeze- 
out characteristics of the hadron phase, making it possible 
for a reliable extraction of the QGP viscosity from exper- 
imental data. 

The default switching temperature T^w to switch hy- 
drodynamics to the hadron cascade simulation is 165 MeV, 
which is chosen from the chemical freeze-out tempera- 
ture measured at RHIC [37] and is approximately close to 
the QCD phase transition temperature from Lattice QCD 
simulations [38] . This is almost the highest temperature 
to implement UrQMD without partonic degrees of freedom. 
It is also the lowest possible temperature for hydrody- 
namic description without introducing additional parame- 
ters for hadronic viscous effects (including both viscosities 
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and relaxations times) and sequential chemical freeze-out 
(which is realized in hydrodynamics by partial chemical 
equilibrium EOS with effective chemical potentials for dif- 
ferent hadron resonances^]- With such default setting for 
VISHNU, the evolution of the system is described as a hy- 
drodynamic expansion of the viscous QGP fluid with a 
chemical equilibrium EOS, followed by microscopic evolu- 
tion of the hadron resonance gas through UrQMD, in which 
the sequential chemical and thermal freeze-out for varies 
hadron species arc realized through the elastic, in-clastic 
and semi-clastic collision rates in the Boltzmann equa- 
tions. 

In Ref. [18]. we investigated whether the microscopic 
hadron cascade approach can be replaced by the macro- 
scopic hydrodynamic approach (with temperature depen- 
dent specific shear viscosity ri/s(T) and partially chem- 
ically equilibrated EoS s95p-PCE as input) by varying 
the switching temperature Ts^j in VISHNU simulations. We 
found that with a constant rj/s as input, the elliptic flow 
shows a strong Tguj-dcpcndencc. After extracting a tem- 
perature dependent effective hadronic shear viscosity (jj/ sY 
from integrated V2 data (with an assumption of short 
relaxation time Tjr = 61] /(sT) from kinetic theory), we 
found that pure viscous hydrodynamics with {rj/ sY^^ {T) 
could nicely fit the pT-spectra and differential elliptic flow 
V2{pt) for identified particles calculated from VISHNU. How- 
ever, the extracted effective hadronic shear viscosity de- 
pends strongly on the pre-hydrodynamic history, particu- 
lary the chosen value of the QGP shear viscosity due to the 
very possibly large relaxation time of the hadronic mat- 
ter. It therefore does not represent the intrinsic transport 
properties of the hadronic matter, but a parameter that 
reflects some memories of the QGP transport properties. 
An extraction of both the shear viscosity and relaxation 
time from UrQMD demands a huge amount of computing 
resources, which is beyond our current scope of investiga- 
tion. Therefore, there exists no switching window below 
Tch in VISHNU, where viscous hydrodynamics can replace 
the hadron cascade [18| . To maintain the predictive power 
of VISHNU, Tsw is to suggested to set at 165 MeV, with 
which the other left free parameters can be fixed from 
experimental data. 



path and instantaneous thermalization, the largest collec- 
tive response is expected. The shear viscosity reduces the 
development of collective flow, leading to a suppression of 
the elliptic and triangular flow as observed by different 

groups [iiaiiniiniiiiiiiiiiisniiiniiii]. 

Viscous hydrodynamics is a useful tool to study the 
viscous effects on the QGP fireball evolution and final 
observables. During the past years, several groups have 
independently developed {2+l)-d llllini[IIl[I21[Il[Il] and 
(3+l)-d 16, 231124] viscous hydrodynamic codes with/without 
longitudinal boost invariance for relativistic heavy-ion col- 
lisions at RHIC and LHC energies. Past research showed 
that shear viscosity decelerates the longitudinal expan- 
sion, but accelerates the transverse expansion, leading to a 
shorter QGP lifetime, more radial flow and flatter hadron 
PT-spectra [TU]. More importantly, it was found that the 
elliptic V2 pi lTUllTTlfT^fT^IH] are very sensitive to the shear 
viscosity. Even the conjectured lower bound from the AdS/CFT 
correspondence rj/s = l/Air leads to a large suppression of 
vo- Thus one can extract the QGP shear viscosity from ex- 
pS^mentsd data by a systematic fltting of V2 as a function 
of collision energy, centrality, system size and etc. 

The flrst attempt to extract the QGP viscosity from 
the elliptic flow data, using 2-|-l-d viscous hydrodynam- 
ics, was done by Luzum and Romatschke around 2008 0. 
They implemented two initial conditions from optical Glauber 
and KLN models and found that the ~ 30% uncertainties 
in initial eccentricity lead to ^ 30% uncertainties for the 
elliptic flow from viscous hydrodynamics with the same 
r]/s, which then translate into ^ 100% uncertainties for 
the extracted value of the QGP shear viscosity. Two ef- 
fects that are neglected in this work are the off-equilibrium 
kinetics (or so-called highly viscous hadronic effects) 
|43] and the partially chemically equilibrated nature [44j 
H5] of the late hadronic evolution, which work against 
each other on influencing V2 and may cancel to some ex- 
tend. Furthermore, initial state fluctuations are neglected 
in their calculations and the effects from bulk viscosity 
was unclear around that time. After making generous es- 
timations for all these uncertainties, it appears that the 
averaged speciflc QGP shear viscosity, cannot exceed the 
following conservative upper limit [9l| 



3 QGP viscosity from elliptic flow data: 
-the early attempt 

Elliptic flow and higher order flow coefficients are impor- 
tant experimental observables for the bulk matter created 
in relativistic heavy ion collisions. In the language of hy- 
drodynamics, pressure gradients convert the initial defor- 
mations and inhomogeneities of the fireball into fluid mo- 
mentum anisotropics at different orders, which translates 
into the asymmetry of particle production as described by 
the flow coefficients. The shear viscosity controls the ef- 
ficiency of this hydrodynamic conversion. In the limit of 
zero shear viscosity, which corresponds to zero mean free 



^ The bulk viscous effects near the phase transition are ne- 
glected here for simplicity. 
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A brief note on triangular flow: 

Recently, several groups extend single-shot hydrody- 
namic simulations to event- by- event ones, making it pos- 
sible to investigate initial state fluctuations and higher or- 
der flow coefficients p?til6 , 35 . 46„47iii5] . It was found that 
triangular flow V3 are also sensitive to the QGP shear vis- 
cosity as the elliptic flow V2 ■ While the initial eccentricities 
£2 differ by 0(20%) between MC-KLN and MC-Glauber 
models, the triangular deformation £3 are almost identical. 
As a result, V3 is much less sensitive to these two initial- 
izations compared with V2. A systematic and combined 
analysis of V2 and V3 together may reduce the initial con- 
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Fig. 1. (Color online) eccentricity-scaled elliptic flow as a function of final multiplicity per area [19] 



ditions ambiguities and give an even accurate extracted 
value of the QGP shear viscosity 0. 



4 QGP viscosity at RHIC and LHC energies: 
-the current status 



sensitive to the form of the non-equilibrium distribution 
function (5/ and the inputting bulk viscosity. In contrast, 
such sensitivity is greatly reduced for the integrated v-i 
(with spectra as a weighted function for the integra- 
tion). We thus proposed to use the integrated for all 
charged hadrons to extract the QGP shear viscosity [T^ 



With the efforts from different groups, the elliptic flow is 
now widely accepted as a key observable to extract the 
QGP shear viscosity. However, it is also significantly af- 
fected by the chemical composition and non-equilibrium 
kinetics of the late hadronic stage as well as the model 
uncertainties of the initialization eccentricity. With the 
newly developed viscous hydrodynamics -I- hadron cas- 
cade hybrid model VISHNU on hand, which more realisti- 
cally describes the hadronic stage and eliminates the re- 
lated hadronic uncertainties, we made an extraction of 
the QGP shear viscosity from the corrected integrated V2 
data at top RHIC energies with emphasis on the remain- 
ing uncertainties related to initialization models and then 
extrapolated our calculation to LHC energies. Below is a 
brief summary of our recent results: 



4.1 QGP viscosity from RHIC integrated vi data 

The hydrodynamic pressure gradients translate the initial 
fireball deformation Ex into fluid momentum anisotropy 
Ep. Meanwhile, the shear viscosity suppresses the devel- 
opment of Ep during the fireball evolution through the 
anisotropic shear stress forces. The experimental observ- 
ables that most directly related to Ep is the integrated 
elliptic flow v'^ for all charged hadrons. While its dis- 
tribution to the differential V2{pt) for identical particles 
strongly depends on the chemical composition and radial 
flow of the hadronic matter, which constantly evolve dur- 
ing the fireball evolution even when Sp reaches saturation. 
Furthermore, V2{pt) at higher pT region {pT > 1 GeV) is 



Currently, VISHNU only concentrate on investigating V2, 
since requires much more computing resources due to the 
event-by-event simulations. 



In Ref. [12], we found that the theoretical v'^ je^ curves 
as a function of multiplicity density per overlap area dNch/(dyS) 
arc approximately universal, which are not very sensitive 
to the initialization models and show clear separations be- 
tween curves as the QGP specific shear viscosity increased 
by l/47r. Furthermore, pre-equilibrium flow and bulk vis- 
cosity only slightly affect these theoretical curves, which 
are at or blow the order of 10%. It thus preferable to ex- 
tract the QGP viscosity from a comparison between the 
theoretical and experimental v'^ je^ — dNch/ (dyS) curves, 
as shown in Fig. 1. the solid and dashed lines with sym- 
bols are the VISHNU results with different {ri/s)QGP as 
input. Left and right panels correspond to two different 
event-averaged, smooth initial conditions, which are ob- 
tained through averaging a large number of fluctuating 
initial entropy densities (given by MC-Glauber or MC- 
KLN models) by aligning the participant plane for each 
event. The experimental flow measurements are generi- 
cally contaminated by non-flow and fluctuation effects to 
some extends, and are not suitable for direct compari- 
son with these theoretical event-averaged V2 ■ Here we use 
the corrected elliptic flow data {v^'^) in the participant 
plane that removes non-flow and fluctuation effects, giv- 
ing an almost universal curve for different corrected flow 
data from most central collision to most peripheral col- 
lisions [in]. While {v"^^) and dNch/dy [SD] are from ex- 
perimental measurements, the theoretical inputs e and S 
are calculated from the MC-Glauber and MC-KLN mod- 
els. This leads to the differences in magnitude and slight 
changes in slope for the two experimental curves shown 
in Fig.l left and right. As a result, the extracted value 
of {r]/s)QGP from these two panels changes by a factor 
of 2 mainly due to the different Ex calculated from MC- 
KLN and MC-Glauber models. Recent event-by-event vis- 
cous hydrodynamic simulations showed that the triangu- 
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lar flow 7J3 is also sensitive to the QGP shear viscosity [TBI 
I15j which suggests that, in the near future, the combined 
analysis of V2 and from VISHNU could yield an even more 
precise value of {f]/ s)qgp and may reduce the ambiguities 
in initial conditions from the hydrodynamic side. At this 
moment, we take the current uncertainties from the ini- 
tial conditions and conclude from Fig.l that the averaged 
specific shear viscosity of the QGP created at top RHIC 
energies is [15] : 



Att s 



< 2.5 X -5- 
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4.2 Estimation of other effects 

In this section, we briefly estimate the residual effects on 
the extracted value of the QGP shear viscosity that are 
not directly included in the VISHNU calculations shown in 
Sec 4.1 and Ref. [19]. 

Bulk viscosity: 

Bulk viscosity also suppresses elliptic flow V2 like shear 
viscosity |13[l29j . Whereas shear viscosity directly sup- 
presses the development of flow anisotropics, bulk viscos- 
ity suppresses the development of radial flow, which in- 
directly influences the elliptic flow through changing the 
slope of the pr spectra and the lifetime of the QGP fire- 
ball. It is generally believed that the bulk viscosity to 
entropy density ratio C/s reaches a peak near the QCD 
phase transition Tc, while the value of this peak is murky, 
for which weakly coupled QCD [51], strongly coupled N=4 
SYM theory [55] and lattice QCD [53] gives dramatically 
different predictions. In Ref [IH], we found that the tem- 
perature dependent bulk relaxation time that describes 
the critical slowing down near the QCD phase transition, 
greatly offsets the effects from the strong growth of C/s 
near Tc, greatly reducing the bulk viscous suppression of 
V2 even when {C/s)max is very large. For simplicity, we 
neglected bulk viscosity in the VISHNU calculations and 
predicted that the extracted values of the QGP specific 
shear viscosity will not be largely contaminated by the 
bulk viscosity due to critical slowing down near Tc- Al- 
though the inclusion of bulk viscous effects could reduce 
the extracted {ii/s)qgp, the total effects will be less than 
20% [29]. 

Event-by-event vs single shot calculation: 

Current VISHNU calculations employed single-shot hy- 
drodynamics with event-averaged, smooth initial condi- 
tions (from MC-Glauber or MC-KLN model) followed by 
thousands of UrQMD simulations to obtain enough statis- 
tics for the spectra and V2 ■ This is a computationally effi- 
cient way to investigate the fluctuation effects on elliptic 
flow, and reduces the computing time by a factor of more 
than 20 when compared with the event-by-event simula- 
tions. However, it also raised the question on how much 
uncertainties it brings to the extracted (rj/s)QGP due to 
the lack of real e-b-e simulations. A detailed comparison 
between e-b-e and single shot simulations from the pure 



hydrodynamics showed that e-b-e simulations reduce V2 
by 0(< 10%) for the same initial eccentricity [15]. This 
indicates that future e-b-e VISHNU simulations will reduce 
the extracted value of {rj/ s)qqp from the elliptic flow data 
by 0(< 30%). 

Initial flow: 



The radial and elliptic flow may develop during the 
early evolution of classical gluon fleld and in the pre- 
equilibrium partonic stage before thermalization. |45ll54j . 
In Ref [T^, we studied the initial flow effects by tuning 
the initial starting time of hydrodynamics with the con- 
straint from fltting the experimental px spectra and flnal 
multiplicities. It turns out that the integrated V2 is maxi- 
mally increased by 0(10%) by the initial flow. This trans- 
lates into an increase of extracted value of {"q/ s)qgp by 
0(< 30%) from this effects alone. 

The form of the viscous correction 5f: 



The form of the viscous correction 5f to the equilib- 
rium distribution is an assumption in viscous hydrody- 
namic calculations. Although it does not directly influ- 
ence the evolution of the QGP flreball, it affects the fl- 
nal pt spectra and differential V2{pt) during the freeze- 
out procedure through the modified Cooper- Frye formula, 
which also infiuences the particle profiles propagated into 
the succeeding UrQMD simulation in the VISHNU hybrid ap- 
proach. Ref. }55j showed that the quadratic and linear 
ansatz of Sf with different assumptions for the relaxation 
time could lead to an obvious difference in the V2{pt) for 
Pt > i GeV. However, the effects at low pr arc small, 
and thus only slightly influence the integrated V2 within 
an order of 5%. Therefore, we choose the integrated V2 as 
a preferable obserables to extract the QGP shear viscos- 

ity m- 

Other initialization models: 



The main uncertainties for the extracted {ri/s)QGP 
are from the undetermined initial conditions from MC- 
Glauber and MC-KLN models used in current VISHNU 
calculations. The fluctuating sources for both models are 
from the fluctuating position distributions of nucleons in 
the colliding nuclei. In Ref. [55] and [57], the additional 
quantum fluctuations for color changes are investigated 
under the framework of Color Glass Condensate(CGC). 
Combing the Classical Yang-Mill's approach for Glasma 
field with the the impact parameter dependent saturation 
model (called IP-Glama model) , Ref. [55] gave an modified 
initial eccentricity £2 and £3 that mostly lie between the 
ones from MC-Glauber and MC-KLN (except for the most 
central and most peripheral region jf). In Ref. [53] , the fluc- 



Ref. [57] investigated the transverse correlations for energy 
density fluctuations within the framework of Color Glass Con- 
densate, but did not further calculate the initial eccentricity, 
since this needs to construct a new Monte-Carlo initialization 
generator to produce the initial energy density profiles with 
correlated fluctuations, rather than un-correlated ones from 
commonly used MC-initialization models. This is still under 
investigation [58) . 
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Fig. 2. (Color online) Transverse momentum spectra of pions and protons for 200 A GcV Au+Au Collisions at different 
centralities [2D]. 



tuations of the initial gluon production are investigated 
accounting to a negative binomial distribution virithin the 
kx factorization approach of CGC. It gives an initial ec- 
centricity £2 that is very close to traditional MC-KLN one 
with only geometry fluctuations, but obvious larger £3 — £5 
than the traditional MC-KLN ones. Considering that our 
current (j]/ s)qgp is extracted from the elliptic flow data 
driven by £2, the further implementation of other initial- 
izations from Ref [SB] and Ref [53] will not change the 
current error bound of {rj/ s)qqp. As discussed in Sec. 3, 
further investigations of elliptic, triangular flow and higher 
order flow harmonics together may help us to distinguish 
which initialization is preferred by the flow data and may 
give an even accurate value of the extracted [rj/ s)qgp- 



To summarize briefly, we extracted the QGP specific 
shear viscosity from the integrated elliptic flow data us- 
ing VISHNU with MC-Glauber and MC-KLN initiahzations 
as input, and found l/(47r) < {ri/s)QGP < 2.5/(47r). The 
width of this range is dominated by uncertainties of ini- 
tial eccentricities from these two initializations. Small bulk 
viscous effects and proper event- by- event hydrodynamical 
evolution of fluctuating initial conditions may slightly re- 
duce the integrated V2- while early flow may slightly in- 
crease it. Although they should be studied in more quanti- 
tative detail, we expect the total uncertainty band trans- 



lated to the extracted value of QGP shear viscosity may 
slightly shift after cancelations. 



4.3 Pt spectra and differential flow for identified 
particles 

After extracting the QGP shear viscosity {ri/s)QGP from 
the integrated V2 data for all charged hadrons, it is impor- 
tant to ensure that VISHNU with {r]/ s)qgp also nicely de- 
scribes the Pt spectra and differential elliptic flow V2{pt) 
for identified hadrons. This was archived and documented 
in Ref [10]. 

Fig. 2 shows the pr spectra for pions and protons at 
200 A GeV Au+Au Collisions. The experimental data 
from STAR and PHENIX collaborations are compared 
with the VISHNU calculations with MC-Glauber or MC- 
KLN initial conditions and different QGP specific shear 
viscosity as input. From most central to semi-peripheral 
collisions, the theoretical curves are insensitive to initial 
conditions and (??/s)qgp and yield an excellent descrip- 
tion of the experimental data. For most peripheral colli- 
sions, VISHNU with finite QGP shear viscosity give a better 
fit of the experimental data, while the ideal fiuid treatment 
shows a slightly steeper spectra for both pions and protons 
due to the insufficient development of radial flow |20] . 

Fig. 3 shows the differential flow V2{pt) for pions and 
protons from 200 A GeV Au-|-Au collisions at different 
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Fig. 3. (Color online) eccentricity-scaled elliptic flow for pions and protons at 200 A GeV Au+Au Collisions [5D]. Solid symbols 
denote reorganized experimental data V2{2}/ {Spart)^ where the measurements i'2{2} are from the STAR collaboration |61] . 
Solid and dashed lines with open symbols are theoretical V2/e from VISHNU. 



centralities. Experimental data are from the STAR col- 
laboration which is obtained from the 2-particle cumu- 
lant method, which measures W2{2} ~ V (^^2)+'^^ with the 
contribution from event-by-event fluctuations cr^^ (where 
0-2^ = {v2y ~ and non-flow effects 6 ^Mi- The 

theoretical curves are calculated from so-called "one-shot" 
VISHNU with the event-averaged initial conditions by align- 
ing the participant plane which corresponds to a smooth 
initial entropy density profile with an eccentricity Epart 
that approximately equals the event averaged eccentricity 
impart)- (The later is suppose to be the driving force for 
the event averaged elliptic flow v(w|)). Since such VISHNU 
calculations are not the real event-by-event simulations, 
they can not be directly compared with the experimental 
U2{2} data. We then compare the theoretical ratio V2/e 

with the experimental ratio V2{2}/ ^ impart) by assuming 

that experimental W2{2} sa \/(v^ ~ (e"^ t) \/ ^^part) with a 
neglect of the non-flow effects (49| . 

Fig. 3 demonstrates that with the {rj/s)QGP extracted 
extracted from the pT integrated V2 for all charged hadrons 
at 200 A GeV Au-|-Au collisions, VISHNU yields an very 
nice description of the differential elliptic flow V2{pt) for 
the identified hadrons (such as pions and protons) at dif- 
ferent centrality bins [50]. The QGP shear viscosity read 
from Fig. 3, {v/s)qgp - (l/47r) for MC-Glauber and 
iv/s)QGP - (2/47r) for MC-KLN actually hit the low- 
est bound of {ii/s)qqp from Fig. 1 for these two initial 



conditions respectively. The slightly lower {r]/s)QGP here 
is due to non-flow effects in the experimental i'2{2} data, 
which give a positive contribution to the measured elliptic 
flow, leading to a slightly lower value of {r]/s)QGP to fit 
the data. 



4.4 An extrapolation to the LHC energies 

The new measurement for 2.76 A TeV Pb-t-Pb collisions at 
the LHC shows a total charged hadron multiplicity den- 
sity that is about a factor of 2.2 higher than the one for 
200 A TeV Au-hAu colhsions at RHIC ESUM], and indi- 
cates a ~ 30% increase in the initial temperature of the 
QGP fireball. Meanwhile the ALICE collaboration at the 
LHC also discovered a ~ 30% increase in the integrated 
V2 and a similar differential V2{pt) when comparing with 
the one measured by STAR at top RHIC energies |64) . 
This raises the question of how the QGP specific viscos- 
ity changes from RHIC to LHC or if one could extract a 
temperature dependent [t]/ s)qqp{T) from the currently 
available experimental data. 

Fig. 4 shows the comparison of the experimental and 
theoretical integrated and differential V2 at 2.76 A TeV 
Pb+Pb collisions and 200 A GeV Au-hAu collision pi] . 
The experimental V2 data are from STAR and ALICE us- 
ing a 4-particle cummulant method, which is supposed 
to measure V2 in the reaction plane under the assump- 
tion of Gaussian fluctuations. The theoretical curves are 
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Fig. 4. (Color online) Integrated (left) and diflFerential (right) W2 for all charged hadrons at different centralities at 2.76 A TeV 
Pb-(-Pb collisions and 200 A GeV Au+Au collision [21]. 



from VISHNU calculations with the event-averaged MC- 
KLN initial conditions created by aligning the reaction 
plane for each event. One finds a nice fit to the experi- 
mental data with {r}/ s)qcp ^ 0.16 at RHIC energies and 
{v/s)qgp ^ 0.20 — 0.24 at LHC energies. Furthermore 
we also find a nice fit of the V2{pt) for pion and pro- 
tons at 2.76 A TeV Pb+Pb coUisions from VISHNU with 
{v/s)qgp- 0.20 m- 

The theoretical calculation shown in Fig. 4 is done by 
inputting a constant [rj/ s)qgp, which indicates that the 
averaged specific shear viscosity (over the space time evo- 
lution of the QGP phase) slightly increases with collision 
energy. However, this does not necessarily mean that QGP 
fluid is more viscous at the higher temperature region 
reached by LHC since one needs to find one temperature- 
depended {r]/s)QGp{T) that fits both the RHIC and LHC 
data rather than assuming two different constant {ti/ s)qgp 
for RHIC and LHC respectively. A further study in Ref [3T] 
shows that one can not uniquely constrain the forms of 
a temperature dependent (77/5) (T) by fitting the spec- 
tra and elliptic fiow at RHIC and LHC energies. Further- 
more V2 becomes more sensitive to the details of the stress 
tensor initialization at LHC energies [SS]. At the current 
stage, no firm conclusion can be drawn on whether or not 
the QGP fluid turns more viscous at the high tempera- 
tures probed by the LHC. 



5 Summary and Concluding Remarks 

In this article, we review recent results from the newly 
developed hybrid code VISHNU [TH] which combines the 
macroscopic viscous hydrodynamic description for the QGP 
fluid with the microscopic hadron cascade model for the 
subsequent evolution of the hadronic stage. Using VISHNU 
with an EOS that implements recent lattice results, we 
made an extraction of the averaged specific QGP shear vis- 
cosity from the integrated V2 data in 200 A GeV Au+Au 
collisions that removes non-flow and fluctuation effects. 



We found that 1 < 47r(?7/s)QGP < 2.5 for the QGP created 
at RHIC, where the width of this range is mainly domi- 
nated by model uncertainties in the initial conditions |19j . 
Compared to the early extraction based on pure viscous 
hydrodynamics [3], this reduces the previous upper limit 
of [rj/ s)qgp by a factor of 2 due to the greatly improved 
description of the hadronic evolution. The (?7/s)qgp ex- 
tracted from the centrality dependence of the integrated 
V2 of all charged hadrons also provides consistent and nice 
description of the pT-spectra and differential elliptic flow 
V2{pt) for charged hadrons as well as identified pions and 
protons over the entire range of collision centralities in 
200 A GeV Au-f Au collisions [20] • After extrapolating to 
the LHC energies, VISHNU also yields a nice description of 
the integrated and differential V2 for Pb+Pb collisions at 
2.76 A GeV with roughly the same averaged QGP shear 
viscosity extracted at RHIC energies [3TJ[35]. 

All of these past VISHNU results are done with single- 
shot hydrodynamic simulations with smooth initial condi- 
tions averaged over thousands of events with a rotation to 
align the " event plane" [T^[^ . This is a computationally 
efficient way to include the effects of fluctuations on the 
elliptic flow using the computationally demanding hybrid 
model VISHNU. However, a detailed comparison between 
event-by-event and single-shot hydrodynamic simulations 
from VISH2+1 shows 0(10%) deviations for the elliptic 
and triangular flow, especially when the QGP viscosity 
approaches zero [15j . Furthermore, higher order flow har- 
monics va, ■ys and vq cannot be realistically described by 
single-shot hydrodynamics through rotating and averag- 
ing the initial fluctuating profiles [TS] . This indicates that 
more accurate extraction of the QGP viscosity from the 
flow data requires the full event- by- event simulations. E-b- 
e simulations from pure viscous hydrodynamic simulations 
revealed that the shear viscosity suppresses the triangular 
flow more than elliptic flow V2 [T?lll6| . This suggests 
that, in the near future, a combined analysis of V2 and 
W3 using the advanced hybrid model VISHNU could yield 
a more precise extraction of the QGP shear viscosity and 
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strongly reduce the uncertainties from the hydrodynamic 
initial conditions. 
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